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Abstract 


Several hull designs are studied with parametric based volume and area estimates to obtain pre- 
liminary hull forms. The volume and area study includes the effects of technologies which manifest 
themselves in the parametric study through stack length requirements. Subsequently, the hull forms 
are studied using a Reynolds Averaged Navier Stokes analysis coupled with a vortex lattice propeller 
design code. Optimization is done through boundary layer control analysis and through studies on 
the effect of blade loading magnitudes and profiles. Several hull form afterbody shapes are studied to 
improve the powering characteristics of a full stern submarine. The hydrodynamic analysis includes 
the total contribution of the rotor, stator, duct, the body hull form and appendages on the submarine 
drag. Corrections are developed and applied to extrapolate the model scale hydrodynamic results 
to a full scale submarine. The end result gives a speed of the hull form which is used to optimize the 
overall submarine. The basic measure of the overall effect of stern taper is the maximum speed of 
the submarine while keeping the same shaft horse power propulsion plant, the same volume and the 
same area requirements. The hydrodynamic powering efficiency is discussed and evaluated using a 
power coefficient. The resulting designs have marked differences in length but similar displacements 
and similar maximum speeds. The results promise great potential for the full stern submarine. 
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Chapter 1 


Introduction 


1.1 Purpose of Thesis 


The history of submarine design has lead to many improvements. From the earliest 
submarine, the Holland, to the newest attack submarine, advances in submarine 
technology usually proceed in methodical steps. Some revolutionary changes include 
the teardrop hull form of the Albacore and the advent of the nuclear age with the 
Nautilus [4]. Perhaps the next revolutionary change centers around the full stern hull 
form with an integrated propulsor. 

Traditional submarines have their machinery weights located relatively aft to pro- 
vide propulsion. A conflicting requirement to this is the necessity to taper the stern 
to prevent separation near the propulsor. A bridge in these two requirements may 
be the use of the full stern hull form. By maintaining a fuller stern further aft, the 
resulting hull form will gain a considerable amount of arrangeable volume. Thus, it 
is possible to move the machinery aft and, thereby, shorten the overall length of the 
submarine. With a shorter hull form and the same installed shaft horse power, the 
submarine may be able to exceed the speed of a conventional designed tapered stern. 


This thesis examines the issues surrounding the full stern design compromises. 
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1.2 Submarine Design Methods 


Stack length describes the requirement that machinery must lie sequentially along 
the longitudinal axis of the submarine. Parametric based estimates use thumb rules 
to relate some known inputs to the size of the submarine. Some common inputs in- 
clude: crew size, number of days underway, number of torpedos, etc. The parametric 
estimates can also include stack length considerations by requiring the most limiting 
of a specified length of a compartment or a minimum compartment volume. 

One method to study the feasibility of a new submarine design uses parametric 
base estimates of area, volume and stack length. A second method uses parametric 
based estimates of weight. In the end, a submarine design must balance weight and 
volume such that the end product is neutrally buoyant. This study considers both 


volume and weight to use the “best” estimates of each methodology. 
1.3. Full Stern Verses Conventional Taper 


The term “taper” refers to the shape of the afterbody as it narrows down from a 
maximum hull diameter to the diameter necessary for smooth transition into the 
propeller regime. Traditional taper designs reach a tailcone half angle of no more 
than about 20 degrees. By limiting the taper transition to a small amount, fluid’ 
separation from the hull is avoided. 

In a conventional design, the hull taper transitions from the maximum diameter 
to the propeller entrance in a length of 2 to 3 hull diameters. A more unconventional 
design, investigated in this thesis and labeled a “full stern”, has a much larger transi- 
tion rate near the propeller. The full stern narrows from the maximum hull diameter 
to the propeller entrance in a longitudinal distance on the order of one hull diameter. 
Using a form of boundary layer control, fluid separation is prevented by accelerating 
the boundary layer into the annular ring at the entrance to the duct. 


One disadvantage to a full stern submarine is the possible loss of propulsor ef- 
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ficiency due to increased frictional and form drag when flow enters the duct. This 
apparent disadvantage is studied further in this thesis. While the efficiency may de- 
crease, the submarine may actually have a faster maximum speed due to the shorter 


hull form and better length to diameter ratio. 
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Chapter 2 


Design Process 


2.1 Baseline Concept Design 


The baseline submarine, used for arrangements comparisons, is a “conventional” sub- 
marine design. The baseline is balanced using parametric relationships originating in 
a course taught by Capt. (ret.) Harry Jackson entitled “Submarine Design Trends.” 
The majority of the relationships are for submarines with 10,000 shaft horsepower 
(SHP) or more. The established baseline falls below the range of some of these 
parametrics. Consequently, in developing the baseline, comparative naval architec- 
ture trends are used to extend the various relationships to a smaller submarine and 
correctly size those items which are not easily scaled [26]. 

The baseline, shown in Figure 2-1, uses existing submarine technologies and near- 
future technologies that could be applied to the next generation submarine. Devel- 
oped to meet the requirements of reference [11], the baseline has a crew of 40 persons, 
four torpedo tubes with space and weight for 20 torpedoes and a diving depth of 1000 
feet. The baseline has allocated space and weight of 400 ft? and 80 ltons to support 
a reconfigurable mission module. The module size meets the maximum for a variety 
of missions. Of the missions considered, the maximum space is required to support 
one platoon of Special Operating Forces plus their support equipment. The maxi- 
mum weight payload is needed in support of anti-submarine warfare to carry a large 


torpedo payload of 30 torpedoes and handling gear beyond the organic capability of 
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Figure 2-1: Advanced Concept Submarine Baseline 


Characteristics 


Submerged Displacement | 3221 lton 
Speed (kts) ZA 
Installed SH P 5900 


Depth (ft) 

Crew 

Endurance Stores 
Number of Weapons 





Table 2.1: Baseline Characteristics Verses the Goal and Threshold Ship Characteristics 


20 torpedoes. 

The propulsion plant of the baseline submarine is a 5500 SH P Pressurized Water 
Reactor (PWR) with a steam system to turn steam turbine generators and steam 
main engines. Power from the main engines transmits through reduction gears, a 
thrust bearing, coupling, shaft seals, etc., to the propeller. The stack length of all of 
these components is estimated to be 60 feet. The estimated stack length is based on 
scaled interpolation from various sources [8], [12], [14], [25]. Figure 2-2 contains the 
estimated stack lengths for the baseline submarine. Table 2.1 provides a summary of 


the baseline design. 
2.2 Full Stern Verses a Tapered Stern for a Submarine 


A full stern design gives considerable flexibility in the submarine internal and external 


arrangements. The fuller stern provides added buoyancy. This buoyancy allows the 
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Figure 2-2: Baseline submarine engine room stack length based on scaled interpolation to 5500 
Saier 


propulsion equipment to be shifted aft, where the limit may become equipment stack 
length, and thus provide more forward space for command and control functions. 
Recent studies have shown that a full stern reduces the submarine length by 5 to 18 
percent while keeping the same equipment and weights within the submarine envelope 
[24], [30]. Additionally, the shorter submarine has a length to diameter ratio closer 
to the optimal of six to one [4]. 

A typical tapered submarine propulsive coefficient (PC) ranges from 0.7 to 0.9. 
There exist a wide range of estimates from various authors [4], [13], [23]. Thurston 
and Evanbar provide a fundamental discussion on the propulsive efficiency of a body 
of revolution and give insight as to why there exist a wide variation [29]. By assuming 
that the traditional naval architecture formulas for a tapered hull apply to a full stern 
(i.e. ignore the separation), a PC can be used to predict the speed of a full stern 
submarine during early parametric studies. The propulsive coefficient decreases from 
about 0.9 to about 0.8 when a full stern is added to a submarine [23]. The PC, defined 


Dy Ca ci is actually greater than 1.0 for a full stern. Section 4.3.4 elaborates. 





15 








MF arian t 4 


Powe ao Q’ 


(i 
— 
CI 
ca 


Figure 2-3: Profile comparison of the baseline and all variants used in the arrangements study. 


2.3 Comparison of Arrangements 


The baseline and all the variants are balanced with the principle characteristics given 
in Table 2.2. A profile comparison of the baseline and all the variants is shown in 
Figure 2-3. 

Several key features are held constant between the baseline and the variants to 
provide a meaningful comparison. Specifically, the shaft horsepower is maintained at 


5500 SHP and the hull diameter is maintained at 28 ft. The Operations Compart- 
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LBP (ft) ZO 232 199 244 
Sub. Disp. (/tons) p22) 3331 2781 2837 
SHP (hp) 9900 9900 5900 59900 
Be ft) 28 28 28 28 
L/D 9.5 8.3 7.9 (ak Sat 
Volume (ft?) 121,200 125,300 117,600 104,700 106,800 
Wet. Surf. (ft?) 19,200 18,750 17,640 15,140 17,800 
Vol/ Wet Surf 6.67 6.91 
ER length (ft) 48.7 S20 
ER volume (ft?) 27,080 17,080 
Propulsive Coef.* 0.81 0.81 
Estimated Max speed (kts) ile Zee 
Stern Shape Tapered | Full Stern | Full Stern | Full Stern Tapered 
Features PWR PWR PWR PWR PWR 
Elect Dr Elect Dr Elect Dr 
Ext Motor | Ext Motor 
* Note: The PC for all full stern variants is a best estimate value which will allow 
the traditional formulas to give a reasonable estimate of speed. The PC for the other 
variants is a calculated estimate from parametric equations. 





Table 2.2: Principle Characteristics of Each Variant 


ment area and weights, as well as the Reactor Compartment volume, are the same in 
the baseline and all variants. 

When a methodical comparison between full stern and tapered stern submarines is 
carefully done, several difficulties present themselves. Foremost, a comparison based 
on volume alone is inaccurate. In other words, the full stern with equal volume as a 
tapered stern can not be arranged internally with the same PWR mechanical drive 
power plant due to stack length requirements. Accordingly, variant 1 requires a ER 
stack length of 60 feet which results in excess ER volume due to the full stern char- 
acteristic. To design a competitive full stern submarine that is physically realizable, 
technologies to reduce the stack length must be included in the design. Consequently, 
this study considers the hull form (full stern or tapered) and various integral tech- 
nologies. These considerations lead to the development of variants 1 through 4. 

Variant 1 through 3 have the same hull form profile characteristics except for the 
parallel mid-body length. A propulsive coefficient of 0.81 is assumed as the value 


for the full stern submarine. The same PC is used in all full stern variants so that 





the changes in maximum speed represent the variation due to the decreased wetted 
surface area and decreased form drag as the parallel mid-body is reduced. The PC 


BHP 


does not represent the traditional definition of PC = CHP’ Instead, these early 


parametric studies use the PC’ as a number which will allow the traditional formulas 


to give a reasonable estimate of speed. Section 4.3.4 elaborates. 


2.4 Variant Descriptions 


2.4.1 Variant 1 —- PWR and Full Stern 


Variant 1 has the same stack length in the ER as the baseline; however, the volume 
is re-shaped into a fuller stern and bow. The Operations Compartment (OPS) length 
is essentially unchanged. The overall length is reduced by 35 feet, due to the finely 
shaped aft ballast tanks and external items of the baseline being placed in the fuller 
hull. The surface displacement remains about the same. The wetted surface area 
decreases by 2% which reduces the frictional drag. At the same time, the PC decreases 


due to a less efficient propulsor resulting in a decrease in maximum speed of 0.2 knots. 
2.4.2 Variant 2 - PWR, Full Stern and Electric Drive 


Instead of a mechanical drive system as considered in the baseline, this variant uses 
an electric drive system. The advantage is a reduction in stack length from 60 feet to 
48 feet and in required volume by 12% based on extrapolations derived from papers 
by Thomas Dade and Jeffery Dutton [8], [12]. 

Figure 2-4 contains the typical stack length assumed for this variant. Both the 
stack length and volume constraints are reached at essentially the same time; the 
volume constraint is most limiting and drives the engine room length to be 48.7 feet. 
The changes made to this variant shorten the baseline by 47 feet and represent a 


technically feasible engine room layout. 
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Figure 2-4: Electric Drive Typical Stack Length used for Variant 2 


2.4.3 Variant 3 —- PWR, Full Stern, and External Electric Drive 


This variant takes variant 2 and places the drive train exterior to the pressure hull. It 
is envisioned that the propeller is actually a permanent magnet rotor and the mount 
and/or duct is a motor stator [1]. Accordingly, the drive train stack length due to 
shaft line components is completely eliminated. However, the free flood percentage of 
envelope displacement is roughly the same. The external drive motor technology has 
been recently demonstrated and is available for incorporation into a future submarine 
[9]. 

The savings in volume and stack length are roughly estimated using previous pa- 
pers that compare and contrast engine room arrangements such as Dade and Dutton. 
The stack length is estimated to be two-thirds of the original engine room length which 


compares well with the removal of the shaft line components from the stack length in 
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the original engine room. As in variant 2, this variant represents a technically feasible 


engine room layout. 


2.4.4 Variant 4— PWR, Tapered Stern, and External Electric Drive 


This variant answers the question: “If the new technologies are placed in a taper hull, 
how would it compare to the best full stern hull?” The reason that this variant is 
less favorable than the equivalent full stern variant is its lack of buoyancy aft. With 
the high power density afforded external electric drive, the propulsion plant occupies 
relatively little volume. And consequently, the submarine tends to sink by the stern. 
Trim lead can be used to balance the longitudinal trim of the submarine, however, 
in this power dense variant, the lead needs to be located forward of the bow unless 
buoyancy is added aft to support the power plant. To bring the variant to a balance 
condition, the engine room volume is increased until ample buoyancy supports the 
power plant and allows reasonable placement of the lead within the hull envelope. 
Figure 2-5 shows the variation in some submarine attributes as the engine room 
volume is increased to bring the placement of the lead within the hull envelope. The 
increased volume aft increases the submarine displacement and increases the length 
beyond that needed if longitudinal balance is not an issue. Accordingly, the maximum 
speed is less than what could be achieve with a smaller volume. Compared to Variant 
3, this variant’s displacement is 56 ltons greater, length is 45 feet longer and maximum 


speed 0.2 knots lower. 
2.5 Conclusion of Full Stern Arrangements Study 


In general, the stack length limitsa full stern submarine such that all of the advantages 
of a full stern design are not realizable without reducing the traditional stack length of 
a propulsion system. Second, the submarine propulsion system and its power density 
plays a significant role in assessing the benefits of a full stern design. Third, the full 


stern design compliments an electric drive submarine, particularly, if the submarine 
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Figure 2-5: The bow is located at 0 ft, increasing x moves aft. The engine room volume determines 
the LCG of the lead. As engine room volume is increased, the aft buoyancy is increased and the 
LCG of the lead can be slid aft. The amount of lead also increases when engine room volume is 
increased since the submarine must submerge. 


has external electric drive. 

A profile comparison of Variant 3 and 4 can be seen in Figure 2-3. The full stern 
variant is 45 feet shorter and 56 /ton lighter. The 22 knot speed of the two variant’s 
is essentially the same. Each variant has unique features. Variant 3 has a full stern 
with its challenges in propulsor design and external electric drive. Variant 4 has more 
typical challenges of current designs with the added challenges of external electric 
drive. Using parametric estimates for the maximum speed, neither variant clearly 
stands out as the best submarine design. 

More detailed hydrodynamic analysis will be done on Variant 3 and 4 in order to 


more accurately quantify their powering characteristics and to discover some useful 
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quantities by which quantitative comparisons can be made. It should be noted that 
Variants 3 and 4 all have the same technologies within their respective hull envelopes. 
As a result, performance differences represent the effects of the basic hull form dif- 
ferences. The particular power plant technologies used are inconsequential as long as 


the power plant technologies eliminate stack length. 
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Chapter 3 


Lifting Surface / RANS 


Optimization 


3.1 Parametric Studies 


To evaluate the hydrodynamic characteristics, a preliminary design study is necessary. 
The preferred level of detail is a lifting line model of the propeller with a quasi- 
empirical estimate of the viscous effects on the body. One tool available at MIT is 
the Propeller Lifting Line (PLL) [7]. This code approximates the downstream helical 
wake as cylindrical and does not model the highly contracting wake expected in a full 
stern design. This difficulty lead to the Ducted Propeller Lifting Line (DPLL) [27]. 
While this code specifically models the contracting wake, it does not capture the full 
viscous nature of the inflow to the propulsor. Work is underway to improve DPLL 
modeling for a full stern bodies [28]. 

The issues above lead to the use of a Reynolds Averaged Navier Stokes (RANS) 
code to model the viscous nature of the inflow which was coupled with a vortex 
lattice Propeller Blade Design (PBD) code to model the propeller influence. The 
procedure commonly annotated as lifting surface / RANS coupling is discussed is 


several references [18], [2]. 
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3.2 Lifting Surface / RANS Evaluation of Sirenian 


In parallel with the arrangement study, a notional full stern body was developed to 
provide a starting basis for future iterations and improvements. The notional full 
stern body is designated Sirenian. It does not represent an optimal design in any 
respects. The Sirenian provides a starting basis to explore possible changes which 


may improve the powering characteristics. 
3.2.1 Physical Geometry 


The Sirenian has a total body length of 11.9 and a body radius of 1.0. The bow is 
elliptical to z/L = 3.0, parallel midbody from 3.0 to 8.9, then two cubic polynomials 
to close the body. The body offset equations are given in equations 3.1 through 3.4. 
The design Reynolds number for the RANS computations is 268829 based on body 
radius. The duct chord is 1.15 with thickness to chord (t/c) of 0.08 and camber to 
chord (f/c) of -0.02 at a 10 degree angle of attack. The trailing edge is located ar/L = 
11.57381821 and r/R = 0.70211869. The duct thickness and camber distribution are 


shown in Figure 3-1. 


Nose - Ellipse, X = 0.0...3.0 
 _ V2-250~ (1.5 ~ X05)" an 
a 2.25 | 
Midbody - Cylinder, X = 3.0...8.9 
=a) (32) 
Stern - Cubic Polynomial, X = 8.9...11.47937081, H = (X — 8.9) 
R = 0.077216975H® — 0.3243113H? —0.018737882H + 1.0 (3.3) 
Closure - Cubic Polynomial, X = 11.47937081...11.899917 


R = —1.47967903.X° + 51.257933.X? — 592.00706.X + 2279.7425 (3.4) 
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Figure 3-1: Duct thickness and camber distributions. 


The stator has 11 blades with the leading edge tip at x/L = 10.4948 and r/R = 
0.857911. The rotor has 6 blades with its leading edge tip at z/Z = 11.0854 and 
r/R = 0.719517. There is no gap between the blades and the duct. The design J, is 


ira2. 
3.2.2 RANS Grid 


The RANS grid is constructed using INMESH code [5] and is shown in Figure 3-2. 
The RANS code used for this analysis is DT'NSA, the axisymmetric version of the 
David Taylor Navier-Stokes codes developed by Gorski [16]. This code is formulated 
with the full Reynolds Averaged Navier Stokes equations. 

The RANS grid is created with nearly constant spacing of the first cell off of the 
body in the normal directions (y+). In order to capture the viscous sublayer of the 
boundary layer profile, a y+ between 1 and 4 is needed. The y+ parameter is a 
function of Reynolds number, which changes along the length of the body. For all 
bodies used throughout this thesis, a y+ of 1 to 4 is maintained along the entire body 


length to accurately capture the flow dynamics. 
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Figure 3-2: The RANS grid used to model the Sirenian body. 


3.2.3 Rotor / Stator design 


The coupled lifting surface / RANS design method requires iteration on many levels. 


Figure 3-3 shows the many iterations paths to reach a fully converged and balanced 


solution. 


The final product of the Sirenian design is shown in Figure 3-4. The stator and 
rotor blade can be seen through the cut-away view of the duct. The Sirenian design 


models a body of revolution; no appendage drag is considered. 


When performing preliminary design, several items are of recurring importance. 


They are: 


e Swirl. The most efficient propulsor will not necessarily have full tangential 


velocity cancellation in the wake. The additional momentum needed to fully 
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Figure 3-3: General procedure for the Navier Stokes / Propeller Blade Design of a vehicle with a 
propulsor. 
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Figure 3-4: Afterbody of the Sirenian propulsor design. 


cancel the tangential velocity will have an associated viscous drag penalty on 
the blades. Conversely, the presence of a hub vortex is undesirable so swirl 


cancellation near the hub is generally desired. 


e Boundary Layer Ingestion. An efficient design will ingest the entire boundary 
layer into the propulsor. There in, the kinetic energy used to build the boundary 


layer is potentially recovered. 


e Wake Profile Aft of the Body. The axial velocity profile left in the wake of 
the body should be as uniform as possible. Ideally, the body would leave a wake 
that is uniform and equal to the inflow velocity. Thus, it would pass through 


the fluid and then leave the fluid “undisturbed.” 
e Drag and Thrust Balance. At the design advance coefficient (J;), the drag 


28 





of the body and appendages must be exactly overcome by the thrust of the 
propulsor. To achieve this balance, the blade loading on the stator and rotor 


must be adjusted while also maintaining the desired amount of swirl cancellation. 


e Duct Loading Fraction. The duct loading can be varied during the design 
stage to evaluate the overall effect on the propulsor efficiency. Performance of 
accelerating and decelerating nozzle arrangements are discussed in many refer- 
ences such as [21]. Additionally, shifting thrust to the duct allows the blades to 


be less loaded and improve the cavitation characteristics of the propeller. 
Lifting Surface / RANS 


To ensure a converged solution is reached, the mass flow through the propulsor and 
the coefficient of thrust and torque (A‘7, Kg) are tracked during the coupled lifting 
surface / RANS procedure. A sample convergence history for a (2%) change in 


advance coefficient is shown in Figure 3-5. 


Normalized Values 
Normalized Values 





25000 25000 
Iteration Number Iteration Number 


Figure 3-5: History of thrust, torque and mass flow rate, and of body and duct drag following a 2% 
perturbation in advance coefficient. 


As would be expected, the viscous drag on the body and on the duct does not 
vary significantly with the small change in the advance coefficient. Conversely, the 
pressure drag makes significant swings, as much as 30 percent. 


For this design, the boundary layer is fully ingested into the duct as shown in 
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Figure 3-6. Additionally, the tangential velocity 1s shown as a contour plot. The 
inflow velocity is 1.0. The large swath of tangential velocity in the wake is a hub 
vortex caused by non-cancelling circulations. Aft of the body centerline, there is a 
strong hub vortex and an axial velocity less than free stream. These indicate that 
this design could be improved through more iterations of the blade circulation. 

The streamline forward and aft of the duct shows how much of the boundary layer 
is ingested in the propulsor. As discussed by Thurston and Evanbar, the most efficient 
design would have the entire boundary layer ingested to recover the kinetic energy in 
the boundary layer [29]. The wake profile should then be as uniform as possible at a 


constant velocity equal to the inflow velocity. 


|UT: -0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 0.20 | 





Figure 3-6: Boundary layer profiles and the duct streamline are shown for a full stern design. 


In reviewing Figure 3-6, the aftmost boundary layer profile indicates that some 
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improvement in the efficiency of the propulsor is possible. The boundary layer has a 
hump greater than the inflow velocity at about r/R equal to 0.5, this hump causes 
proportionally more kinetic energy to be left in the wake. Kinetic energy is propor- 
tional to V+ V, whereas thrust is proportional to V, thus exceeding the inflow velocity 
tends to waste more energy. Similarly, in Figure 3-6, there is a depression in the wake 
profile near r/R of zero. The depression indicates that by directing more flow aft, 
and incidentally less in swirl, the efficiency could also be improved. 

Improvements to the design are easily identified. However, actually accomplishing 
these improvements is tedious and time-consuming. The improvements are largely 


based on engineering judgment and checked with detailed numerical solutions. 
3.3. Model Scale Powering Characteristics 


Recall this propulsor / body combination does not represent an arrangeable submarine 
as developed in Chapter 2. Rather, this combination is just a representative estimate 
of a submarine. It does not include appendage drag and the frictional drag has not 


been adjusted to a full scale value. The powering results are shown in Table 3.1. 


J Siren 
Cpenre : 
Cie 


PercentOverpowered 
Kr 

10 * Kg 
PowerCoe f ficient 





Table 3.1: Results for the Sirenian model scale design. 
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Chapter 4 


Full Scale Design of 
Sirenian 


4.1 RAWNS Grid 


The RANS grid has been shortened compared to the original Sirenian grid. The main 
reason for this is two fold. First, since comparing propulsor geometries affects the 
afterbody only and having an identical bow does not change the inflow properties, the 
grid size can be reduced to save computational time. Secondly, variations originating 
at the bow due to numerical instabilities could lead to errors in the comparisons. The 
pressure drag tends to behave erratically and does not consistently settle to stable 
value. Figure 4-1 shows the variations of pressure on the body and the duct over 
1000 RANS iteration cycles. The afterbody evaluation of pressure seems consistent. 
However, the evaluation of pressure near the bow and, thus, the evaluated pressure 
drag is unusable. These errors originate from numerical instabilities induced from the 
RANS grid spacing across the zonal junction near the bow. 

First, to solve this difficulty, the body pressure drag is evaluated using a geomet- 
rically generated RANS grid on a similar body. Secondly, to provide an independent 
check, the body pressure drag is evaluated using an Integral Boundary Layer method. 
The results are encouraging. Figure 4-1 shows the comparison. For full scale Sirenian 
and for later variants, the RANS solution solves the afterbody portion of the flow 


only and uses the stable, geometrical gridded, solution as a boundary condition at 
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5 

x/R 
Figure 4-1: The figure on the left shows the body pressure variations over 1000 RANS cycles in 
a stabilized flow field. The figure on the nght shows the resultant body pressure drag integrated 


over the body for identical inflow and bow/midbody geometry. Variations between the methods are 
expected near the propulsor due to different assumed propulsors and stern fullness. 


the upstream edge of the grid. 

Another difficulty presents itself when establishing a RANS grid for the first time. 
The boundary layer tends to develop a discontinuity across the grid junction if the 
grid spacing and grid orthogonality are not matched closely across the grid junction. 
The difficulty can be minimized by matching the grid vertical spacing across the 
junction along the entire boundary axial extent. 

The boundary layer discontinuity can be made a null issue by ensuring the grid 
boundary follows a streamline. This ensures that no fluid needs to cross the boundary. 
The streamline is not known a priori. Consequently, using an open water case stream- 
line then replacing it with the new ducted case streamline provides a mechanism to 
alleviate the boundary layer discontinuity. It should be noted that the streamline 
location changes when the thrust is changed. Thus, to a small extent, the grid zonal 


boundary is iterative with the thrust. 
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4.2 Full Scale Estimates 


In Chapter 3, Sirenian matches drag and thrust on the model scale after solving the 
vortex lattice / RANS coupling. To allow Sirenian to be more representative of a 
real submarine, the effects of model verses full scale are considered. Using the ITTC 
friction line estimate, the model scale and full scale frictional drag are calculated and, 
after correctly non-dimensionalizing the quantities, are used to scale the Sirenian to 


full scale. Equation 4.1 shows the relationship. 
CDsmp = CDsoaee — Cos(Roodel) + Cos (Rsaip) + Cdanowance + CDappendage (4-1) 
The model drag is summed in the calculations as in equation 4.2. 
CDyoaet = CDaprcs + CDaps>s + CDavres + CDav.er8 + CDpp + ODay + CDstator (4-2) 


For different designs, the frictional drag at full scale will vary. However, the changes 
in the actual frictional coefficient are neglibly small. To illustrate by example, the 
assumed full scale condition for all Sirenian frictional coefficients is 199 feet and 22 
knots. A change of 1 knot in speed results in a 0.6 % difference in Cps(Rsnip). 
Similarly, 5 knots results in 2.6 % difference. Recall, the Sirenian is not scaled to 
correspond to the variants of the arrangement study; these will be examined in the 
following chapters. 

A correlation allowance of 0.0005 (non-dimensionalized on wetted surface area) is 
used for all variants. The appendage drag is calculated using Hoener’s relationship as 
discussed in [23]. The appendage drag is expected to add a fixed amount to the Cp 
of 0.022. This additional drag represents about 15% to 20% additional appendage 
drag and is consistent with other references such as Burcher and Rydill [4]. Table 4.1 


summarizes the corrections used to extrapolate the model scale drag to full scale. 
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Cos (Rm oaet) 
Cos (Rsnip) 


C@pallowance 
Cp appendage 
Total Correction 





Table 4.1: Corrections used to extrapolate the model scale Sirenian drag to full scale drag. 


4.3 Powering Characteristics of a Full Scale Sirenian 


The ultimate goal of a design such as this one is to maximize the speed with the 
install power. Assuming no hub, and therefore nearly constant diameter helices in 
the wake, the optimal circulation distribution on a blade set in a duct is constant 
distribution. The upstream blade set induces a circulation which is exactly canceled 
by the downstream blade set. Reference [20] discusses this further. 

For the full stern, the helical wake contracts significantly and has a pronounced 
strong interaction with the rest of the propulsor elements. A constant circulation is 
used as a starting point for this design. Variations from constant are used to explore 


other designs. 
4.3.1 Circulation Distribution 


For all Sirenian designs, the circulation distribution is linear. The slope of the dis- 
tributions vary to provide a range of cases to evaluate. The constant case uses the 
non-dimensional circulation of 0.0110 on the rotor and 0.0110 * (6/11) on the stator. 
This ratio of rotor loading to stator loading provides nearly complete swirl cancella- 
tion. 

The tip loading in excess of the r/R = 0.7 value (0.0110) is used to gage the 
various cases. When varying the circulation, the stator value is adjusted by the same 
tip loading excess times the factor 6/11 to account for the different number of blades. 
Finally, the magnitude of blade loading is adjusted to make the overall thrust and 


drag equal within a few percent. Adjusting the magnitude made the percent of actual 
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excess loading change when the distribution is uniformly scaled to achieve a balanced 


thrust and drag. Figure 4-2 illustrates these ideas. 
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Figure 4-2: The left most figure shows representative non-dimensional circulation distributions. The 
figure on the right denotes how the actual value of non-dimensional circulation at the 0.7 span 
changes to maintain a balanced thrust and drag condition. 


4.3.2 Results of Circulation Distribution Variations 


Table 4.2 summarizes the evaluated cases. The thrust and drag are balanced within 
+ 1.0% for all cases. Note, however, the thrust and the drag both increase as the 
hub is loaded more. Since higher hub loading causes higher velocities close to the 
body, the frictional drag is increased. This is analogous to a higher thrust deduction. 
Additionally, the higher hub loading tends to produce a hub vortex which adds to 
the body drag. The combination of more frictional drag and more hub vortex drag 
is offset by the higher hub loading. Increasing the hub loading, provides more thrust 
with less penalty in the require torque since the radius from the hub is smaller. ‘Thus, 


more hub load provides an overall benefit to the ship’s powering characteristics. 
4.3.3 Conclusion on Blade Loading 


The minimum power for a given case is directly related to the torque coefficient, Ag. 


Since P = 27nQ and using the definitions of A’, Ng and J;, the case with the best 
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[| Cisne | Crenrp | Overpowered | Kr_| 10* Ke 


20%G 
10%G 
5%G 
Constant 
5%L 
10%L 


15%L 
20%L 
30% L 
40%L 
50% L 
60%L 
TO%L 


0.08916 
0.08907 
0.09048 
0.09136 
0.09162 
0.09245 
0.09326 
0.09357 
0.09596 
0.09742 
0.09899 
0.10043 
0.10203 


0.08836 
0.08972 
0.09080 
0.09188 
0.09092 
0.09160 
0.09324 
0.09368 
0.09632 
0.09760 
0.09900 
0.10032 
0.10188 





Table 4.2: The results associated with the non-dimensional Sirenian body / propulsor combination 
indicate higher hub loading provides better powering characteristics. 


speed is the same case that has a minimum value of Ag. The trend clearly shows 
that less tip loading gives a better design when considering speed for a given installed 
horsepower. Some variations in the trend are related to small deviations in the vortex 
lattice / RANS coupling process. Yet, the overall trend is evident and is graphically 
shown in Figure 4-3. 

The trends seems to continue past the values evaluated. The difficulty becomes 
the excessive lift coefficient required to generate the specified hub loading. Using a 
typical lift coefficient upper limit of 0.4 to 0.5, the 60% less tip loading case seems to 
reach the reasonable design limit. Figure 4-4 shows the lift coefficient verses span for 
the 60% less tip loading case. 

Reviewing the axi-symmetric flow field from the RANS output gives the reason 
for the trend toward higher “efhiciency” with higher hub loading. The velocity in the 
wake of the body tends to get filled in by higher velocity fluid. The wake deficient 
velocities are brought closer to the free steam velocity. As the wake is filled in, the 
wake profile approaches a uniform, constant velocity. In light of the discussion in 
section 3.2.3, the trend that higher hub loading gives higher efficiency is reasonable. 

The trend can also be evaluated using the wake deficient analogy. The body 


produces a thick boundary layer upstream of the propulsor. As this boundary layer 
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Figure 4-3: The variation of Cp shows higher hub loading minimizes the power. 


approaches the propulsor, the fluid closest to the body has the most kinetic energy. 
Thus, to recover this energy, the fluid closest to the body needs accelerating the most. 
Hence, higher loading at the hub is needed to recover the energy in the boundary layer. 

Finally, by loading more heavily at the hub, the tip loading is reduced. This has the 
benefit of reducing likelihood of tip cavitation. While this conclusion finds that 60% 
nominal less loading at the tip is a good design, it is not “optimal.” More studies 
are needed to confirm that this trend holds given other constraints. Additionally, 


variations from the linear distribution may yield even better results. 


4.3.4  Propulsive Coefficient 


Now that the propulsor characteristics are calculated, can a propulsive coefficient 


(PC) be calculated? In short to calculate a PC, the unpropelled drag on the towed 
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Figure 4-4: For the 60% less tip loading case, the coefficient of lft on the blade approaches 0.5. 
Separation will likely occur at Cz, higher than about 0.4 to 0.5. 


body must first be estimated. The difficulty here lies in the fact that the unpropelled 
body will exhibit separation near the stern, and present day numerical methods can- 
not accurately model the separated flow dynamics. 

The propulsive coefficient measures the effective horsepower out of the system 
divided by the shaft horsepower put into the system. The effective horsepower is 
towed resistance times the velocity of the ship. The towed resistance, Rr, is related 


to the propeller thrust by a thrust deduction factor [21]. 


EHP 
PC = Fae (4.3) 
eae naa (4.4) 
Rr =(1—t)T (4.5) 
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The towed resistance (i.e. unpropelled resistance) for a full stern is quite large due 
to separation on the after part of the ship when compared to its propelled resistance. 
Separation is prevented by the action of the propulsor. Generally, when the separation 
is removed by the action of the propulsor, the hull drag will decrease. The effect of a 
large unpropelled resistance will make the E'H P larger for a given speed. Accordingly, 
a full stern PC will be higher than a conventionally tapered stern PC’. While the 
PC will be higher, this does not lead to the conclusion that the maximum speed will 
be higher. 

The maximum PC attainable can be larger than 1.0 [29], [21], [27]. Neglecting 
shaft losses (75 = 1.0), the PC for the 60% less tip load case is 1.41(1 —¢t). Using 
equation 4.5, we know that (1 —t) > 1. Thus, the PC will be some value greater 
than 1.41. This does not lead to the conclusion that this design will have a higher 
speed than a design with a lower PC’. The conclusion is that the propulsor in this 
design does a better job of propelling the body verses the body being towed. 

Using a PC for a full stern design is very misleading. ‘The speed, thrust, and torque 
can be calculated directly from the coefficients obtained in a coupled lifting surface 
/ RANS analysis. The best measure of full stern performance is the maximum speed 
attainable for a given installed horsepower. Understandably, however, an “efficiency” 
is desired to use during comparisons. One “efhciency” that measures power is the 
Power Coefficient, Cp. The Power Coefficient is defined in equation 4.6 and is shown 


in Table 3.1 and Figure 4-4 for this design. 


_ 2xnQ 
sPVEAR 


Cp (4.6) 


Each of the variables are defined in the nomenclature section. The Cp represents the 
power input to the propeller divided by a power associated with the ship speed and 


the frontal area of the body. 
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Chapter 5 


Lifting Surface / RANS Evaluation 
of Variant 3 


The Sirenian propulsor attributes are combined with the results from the arrangement 
study to form a physically realizable submarine variant. In order to estimate the 
maximum speed attainable from the installed 5500 SHP, the variant is evaluated for 
Ky and Kg. The trends of Chapter 4 are applied to this full stern variant. A sampling 
of cases are used to evaluate the maximum hub loading without separation. Then, 
the blade loading is adjusted to refine the boundary layer in order to maximize speed. 
Finally, various body stern geometries are used to reduce the drag in order to increase 


speed. 
5.1 Physical Geometry 


The Variant 3 has a total body length of 14.2 and a body radius of 1.0. The bow is 
identical to the bow used for Sirenian. This similarity aids in direct comparisons of 
propulsor changes. The bow is elliptical to z/L = 3.0, parallel midbody from 3.0 to 
11.2, then the same two cubic polynomials used in the Sirenian design close the body. 
Other afterbody shapes are discussed in section 5.4. The design Reynolds number 
for the RANS computations is 268829 based on body radius. The duct chord is 1.15 
with t/c of 0.08 and f/c of -0.02 at a 10 degree angle of attack. The duct has been 


moved inward to capture less of the free stream and place the leading edge at the 
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[] Cieme | Creare | Overpowered | Kr | 10* Kg 


20%G chez. 0.12192 


Constant 0.12568 


20%L 0.13076 
40%L 137 0.13676 
60%L 0.14385 





Table 5.1: The results associated with the non-dimensional Variant 3 indicate that the nominal 20% 
less tip loading case has the best powering characteristics. 


| eT healt) eee MazSpeed(knt) 


20%G 98142. 340601. 
Constant 102380. : 338580. 


20%L 107076. 337698. 
40% L 111648. : 338213. 
60%L 116319. 339821. 





Table 5.2: The results assume a submarine length of 199 feet and 5500 SH P. 


“knee” of the boundary layer. 

The stator has 11 blades with the leading edge tip at #/Z = 12.7948 and r/R = 
0.740817. The rotor has six blades with its leading edge tip at ¢/LZ = 13.3854 and 
r/R = 0.602437. There is no gap between the blades and the duct. The design J, is 
1.52. The selected RANS grid is similar to the RANS grid used to calculate the Full 


Scale Sirenian. The grid starts at x=8.0 and extends to x=80. 
5.2 Results of Circulation Distribution Variations 


Table 5.2 assumes a physical submarine that matches Variant 3 from the Chapter 2 
arrangements study. Since the length to diameter ratio of Variant 3 is already selected 
to match the Chapter 2 study, the only additional assumptions are a length of 199 ft 
and an installed power plant of 5500 SH P. Using a linear blade loading distribution, 
a gain of 0.2 knots on the maximum speed is possible when the tip has 20% less 
nominal tip loading. The design process undertaken uses the propeller torque and 
RPM as variables, not fixed at the onset. This is consistent with the use of electric 


drive with solid state speed controllers as used in the arrangements study of Chapter 


Ze 
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Figure 5-1: The trend shows that a higher hub loading minimizes the power until the loading is so 
high that other influences such as frictional drag, hub vortex, etc. offset the gains. 
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Figure 5-2: For the 20% less tip loading case, the coefficient of lift on the blade approaches the 
highest that can be reasonable obtained. As the C, approaches 0.5, separation will likely occur at 
the rotor hub. 
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The local optimal of the 20% less tip loading case has the lowest power coefficient 
and thus the highest speed as shown in Figure 5-1. Additionally, this case has the 
highest rotor lift coefficient that can be reasonably obtained. Figure 5-2 shows the 


lift coefficient for this case. 
5.3 Blade Loading at Hub 


Thus far, more hub loading gives better powering characteristics. However, when the 
blade loading causes the blade section coefficient of lift to exceed a reasonable value, 
the cavitation and flow separation on the blade is likely. To capitalize on these trends, 
the blade loading is adjusted from a linear distribution to a distribution which loads 
the blade hub more without exceeding a lift coefficient of about 0.45. 

The revised blade loading uses the 20% less tip loading case and then, on the 
blade sections less than about 30% span, the blade loading is reduced. Finally, the 
blade loading is adjusted so that a reasonable coefficient of lift is reached on the lower 
30% of the blade. The blade loading and the resulting lift coefficient are shown in 


Figure 5-3. Table 5.3 gives the powering characteristics. 
5.4 Body Stern Geometries 


Figure 5-4 shows the integrated drag due to pressure and viscous influences on the 
two-cubic polynomial Sirenian afterbody. Considerable drag occurs on the body near 
the propulsor entrance way. The body and the duct have a strong interaction on each 
other’s drag components. The interaction of the duct on the body plays a role in the 
total drag on the body. In this case, the body drag is considerable while the duct 
actually provides a net thrust. 

The effects of body profiles are studied to better understand the sharp rise in drag 
near the propulsor. The boundary layer is highly sheered as it is accelerated towards 


and through the propulsor. Due to the sheered flow near the propulsor, an increase 
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Figure 5-3: The best results load the hub so that the lower one-third of the blade has a lift coeffi- 
cient that approaches the maximum attainable without blade flow separation. Other blade loading 
distributions which build on these trends could provide even better results. 


in the viscous drag is expected. A corresponding increase in the pressure drag is also 
expected, although the anticipated magnitude and shape are less clear. 

To mitigate the increase in drag near the propulsor, a smoother transition is needed 
as the flow is redirected from the parallel midbody into the propulsor. The chosen 
functions, y = z/(1 +2) and y = tanh(z), are used to provide a general shape. 
To cast the functions into the shape of a full stern, the functions are multiplied 
by constants so that their shape matches the original Sirenian shape at the stern. 
Finally the functions are joined with a third-order polynomial to close the stern. For 
all geometries, the transitions maintain a continuous first and second derivative to 
alleviate pressure discontinuities in the flow. Figure 5-5 shows the three body profiles 
and their derivatives. 

From the cases studied, the geometry which has the smoothest transition, 1.e. 
continuous higher order derivatives, is also the geometry which has the best pow- 


ering characteristics. The fractional case, y = z/(1 +2), has 21% less drag than 
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Figure 5-4: The integrated drag due to pressure and viscous influences are shown for the Sirenian 
afterbody. One method to improve the powering characteristics is to reduce the sharp rise in the 
drag at the entrance of the propulsor. 
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Kr 


10* Kg 
PowerCoef ficient 
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ShaftRPM 
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Table 5.3: The smoother the body profiles, the better the powering characteristics. 


the cubic polynomial geometry of the Sirenian afterbody. The hyperbolic tangent 
case, y = tanh(zx), has 22% less drag than the Sirenian afterbody. The later case 
has the smoothest transition from the parallel midbody to the propulsor. Table 5.3 


summarizes the powering characteristics. 
5.5 New Body Geometries 


The general trends need to be identified to establish the characteristics of a optimum 
powering full stern body. The techniques used throughout this thesis can be used to 


establish the trends. By studying the flow details around a series of bodies, the trends 
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(a) The body profiles are shown with an exag- 
gerated ordinate scale. The bold line and open 
square symbols show the original body profile. 
The fractional body is shown with right trian- 
gles. The hyperbolic tangent body is shown 
with a dotted line and circles. Note the close 
similarity in the body profiles. 
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(c) This version uses the curve y = £/(1+ 2) 
to provide the curvature near the propulsor. 
A third-order polynomial is used from 13.9 to 
14.2 to close the body. The first and second 
derivatives are smoother than the original Sire- 
nian body. 
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(b) The Sirenian body profile and the corre- 
sponding first and second derivatives are shown 
above. The Sirenian uses two third-order poly- 
nomials to represent the body in the propulsor 
region. 


Body Offsets and Derivatives 
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10 11 
x/R 
(d) This version uses a y = tanh(x) curve to 
provide the curvature near the propulsor. A 
third-order polynomial is used from 13.9 to 
14.2 to close the body. The first and second 
derivatives are smoother than the other bodies. 
This version has an infinite number of contin- 
uous derivatives in the region on z/L = 11, 
which is where a great reduction is achieved in 
both body pressure and viscous drag. 


Figure 5-5: Body profiles and derivatives 
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can be developed and applied to a parent body to establish an optimum design. 

Future studies can use the results contained herein to establish a parent hull form. 
The new parent hull form could serve as a central point for the variations of each 
characteristic. As long as the tested characteristics are close to the parent hull form 
characteristic, the results should extend to the case where several characteristics are 
simultaneously changed to their respective optimal value. While this idea does not 
optimize the whole body at once, it does give a reasonable approach to an extremely 
broad problem. 

In order to ensure continuity in higher derivatives to minimize drag, the parent hull 
form should be constructed out of a single mathematical function. One possibility 
is the use of B-splines. B-splines describe the body with much fewer control points. 
The order of the B-spline curve is determined by the number of constraints placed on 
the B-spline. 

A second possibility is the use of a power series of the form, y? = a,z+a9z7+...+ 
a,x”, where x is the non-dimensional abscissa and y is the non-dimensional ordinate. 
The arbitrary constants a), a2, etc. for each form are determined by the geometrical 
properties of the form. The number of constants is determined by the number of 
constraints on the body. 

The constraints include the body properties such as y(L) = 0 and the naval ar- 
chitecture characteristics such as the aft prismatic coefficient C'p,. The following list 


shows possible constraints for a full stern body series. 


e Assuming the power series describes only the curved afterbody, continuity of 
derivatives should be maintained by specifying position and derivatives at the 
parallel midbody and afterbody junction. As a minimum, the slope and curva- 
ture should be matched. However, matching higher derivatives is expected to 
improve powering. An alternative method describes both the parallel midbody 


and the afterbody with the power series. This method still requires specify- 
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ing derivatives to ensure the parallel midbody is maintained up to a specified 


location. 


The radius of curvature can be specified at the three primary locations on a full 
stern body. Moving aft along a typical full stern afterbody, the curvature will 
first have a large negative maximum at the shoulder, a large positive maximum 
somewhere in the middle and, finally, a large negative maximum as the body 
closes. Both the (z,y) location and the magnitude of maximum curvature can 
be specified for each location. The radius of curvature is given by: 


eee lal 


p(x) d? y(z) 
dx? 
The location could be specified by y(z,) = y,; where zx, is determined from 


dp(x1) oa 0 
df 


The prismatic coefficient can be specified in the power series representation. The 
COMstraint Is Of the 1oLml. Cp, — REL Si, y(x)*dx where L, is the length of the 


afterbody. 


The length to diameter ratio can also be specified by the location of the tail 


cone. 


The duct characteristics must be somehow integrated in the body series study. 
The duct characteristics which could be varied include: angle of attack, thickness 
and thickness distribution, camber and camber distribution, location of duct with 
respect to the body, location of the duct with respect to the rotor and the radius 


of the duct inner surface. 


Since the full stern body must be evaluated while under power to capture the 
dynamics of the boundary layer, characteristics of the propulsor blades must 
also be included in a body series. These characteristics include: blade loading 


distribution and division of loading between the stator and rotor. 
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Chapter 6 


Lifting Surface / RANS Evaluation 
of Variant 4 


The next comparison evaluates the arrangements variant 4. This case represents a 
more typical submarine afterbody shape than the full stern Variant 3. The typical case 
uses the same elliptical bow as all other variants and uses the Huang 1 afterbody [17]. 
The overall length is achieved by varying the amount of parallel midbody to match the 
Variant 4 length to diameter requirement. Both a non-ducted and a ducted version 


are evaluated. 


6.1 Non-Ducted Tapered Body Evaluation 


6.1.1 Physical Geometry 


The Huang 1 afterbody is used to represent typical tapered stern submarine geometry, 
as well as typical propeller placement and characteristics. Since this case represents a 
typical case, the given geometry and placement are taken as optimum given realistic 
constraints on the design. 

Figure 6-1 shows a notional sketch of the geometry. Consistent with reference [17], 
the design has seven blades with the leading edge tip at x/L = 17.1 anda r/R at the 
tip of 0.545. No duct nor pre-swirl vanes are present. The design J, is 1.25. 

Using similar methods to those used in section 4.2, the full scale corrections were 


developed and are shown in Table 6.1. 
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Figure 6-1: Notional tapered stern non-ducted design 





















Cos (Ra odet) -0.09944 
Cos (RsurpP) 0.04481 
C Daliowance 0.01446 
C'Dappendage 0.0220 


Total Correction } -0.01817 


Table 6.1: Corrections used to extrapolate the tapered stern results to full scale drag 


ol 





[Tapered with No Duct 


CDarane 

CTsutp 
PercentOverpowered 
Kr 

10 * Ka 


PowerCoe f ficient 
Thrust (lbf) 
ShaftRPM 
Torque(ft —lbf) 
MazSpeed({knots) 





Table 6.2: Results for the open water propeller 


6.1.2 Powering Characteristics 


The circulation distribution used for the open water case is assumed to be an optimum 
distribution based on previous work [17], [20]. The powering results are shown in 
Table 6.2. By using the same technologies for this variant as the full stern variant, 
neither the torque nor the shaft speed have limits imposed. Consequently, they are 


free to float to the values that best optimized the propulsor design. 


6.2 Ducted Tapered Body Evaluation 


6.2.1 Physical Geometry 


The Variant 4 ducted design uses as much in common with the non-ducted Variant 
4 and with the full stern Variant 3, as possible. As such, the duct chord is 1.15 with 
t/c of 0.08. However, due to massive separation on the duct outward side when the 
camber and angle of attack matched the original values, the f/c is 0.01 with a 6 
degree angle of attack. To reduce the separation on the duct, the inner surface of 
the duct is aligned with the streamline in the open-water propelled case that passed 
through the location of the propeller tip. The rotor tip is placed at the same axial 
location on the body as in the non-ducted case. The reduced camber and reduced 
angle of attack will cause the duct to generate less lift in the direction of forward 


tomas. 
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The radial placement of the duct inner surface is place using the ideas presented 
in reference [15]. The placement maintains the same thrust between the open water 
propeller and the ducted propulsor. The r/f of the rotor is 0.470 and places the 
ducted inside the “knee” of the boundary layer. This judgment observation is clouded 
by the fact that the boundary layer diffuses over the tapered stern, and thus the radial 

location of the “knee” changes with axial position. 

The stator has 11 blades with the leading edge tip at x/L = 16.5094 and r/R = 
0.5267894. The rotor has 6 blades with its leading edge tip at 2/Z = 17.1 and 
r/R = 0.470. There is no gap between the blades and the duct. The design J, is 
1.52. The number of stator blades, the number of rotor blades and design J, are 


consistent with the full stern ducted design. Figure 6-2 shows the notional geometry. 
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Figure 6-2: Notional tapered stern ducted design 


The RANS grid models the same flow field area as the full stern grid. The grid 
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[~*~ Papered with Duct 


CDsute 0.09018 
Crone 0.08994 
PercentOverpowered -0.3 

Kr 0.3694 
10* Kg 0.8172 


PowerCoe f ficient 0.0822 
Thrust (lbf) 84490 
ShaftRPM 117.4 
Torque 245973 
MazSpeed(knots) Zee 





Table 6.3: Results for the ducted, tapered stern submarine. 


starts at x/L = 8.0 and extends to z/Z = 80.0. The number of cells is maintained 
similar to the number of cells in the full stern grid. Therefore, the grid resolution is 
very nearly identical to the full stern grid. As in the non-ducted case, the full scale 


corrections are shown in Table 6.1. 
6.2.2 Powering Characteristics 


For a counter rotating set of propeller blades in a non-contracting flow field such as 
behind a gradually tapered afterbody, the optimum circulation distribution is very 
nearly constant loading [17], [20]. The powering results are shown in Table 6.3. Just 
as the other cases, neither the torque nor the shaft speed have limits imposed. 
When compared to the open water version, the ducted version of Variant 4 demon- 
strates the advantage of a ducted design. With a duct present, the Variant 4 gains 0.9 
knots in maximum speed. While there are literally hundreds of variables to adjust to 
obtain the optimum design, this simple excursion shows clear advantages of a ducted 
design over a open water design. These ideas are not new, they have been discussed 


in countless papers and reports including [20], [29], [10], [3]. 
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Chapter 7 


Conclusion and Recommendations 


7.1 Summary of Results and Conclusion 


The design process followed in Chapter 2 yields an honest comparison of a full stern 
and a tapered stern hull envelope. The parametric arrangements, while not exact, 
provide a good comparison of two similar versions of the same submarine. The pow- 
ering characteristics developed in Chapters 3 through 6 provide a comparison of the 
propulsor performance which corresponds to the submarine hull envelopes developed 
in Chapter 2. 

The ultimate goal of this thesis is to answer the question of whether a full stern is 
a better hull envelope than a tapered stern. Many assumptions are made throughout 
this study to obtain quantifiable results. Whenever possible, the assumptions are 
forward looking while still realistic. The hope is to include the effects of interior 
arrangements on the size of the envelope. Then, using the correctly sized envelope, 
determine the maximum speed attainable with the same install shaft horsepower. 

Starting with arrangements and ending with hydrodynamic powering calculations 
yield significant insight on the full stern. The results, while not conclusive, promise 
great potential for the full stern submarine. The two designs have marked differ- 
ences in length but similar displacements and similar maximum speeds. Table 7.1 
summarizes the competing designs. 


The full stern has a faster maximum speed than a tapered stern. The full stern 


29 





speed is 23.5 knots which compares with a tapered stern speed of 23.2 knots. The 
speed difference is only 1.3 % which is well within the precision of the vortex lattice / 
RANS method used to analyze the powering characteristics. However, recall that the 
full stern is not globally optimized. Many more cases can be envisioned which may 
yield even better performance. Thus, the closeness in speed for such a small sampling 
of cases promises great benefits from a full stern submarine. Therefore, additional 
study should be directed towards the full stern submarine to verify this study and 


provide more insight on the benefits and drawbacks of a full stern submarine. 
7.2 Conclusion on Arrangements 


The parametrics used to arrange the submarine have limitations with the issues of 
stack length and with volume to area conversions. Reference [11] does preliminary 
arrangements of a submarine which was initially sized with the parametric arrange- 
ments and gives a conclusion on the precision of the parametrics. The final design 
in reference {11] shows that parametrics yielded at hull envelope of 181 feet and the 
preliminary arrangements yielded 178 feet. Thus, the arrangement parametrics used 
to size the hull envelopes give reasonable size hulls. There was, however, substantial 
differences in the division of volume interior of the hull envelope. 

The arrangements give a hull form envelope on which to base viscous and pressure 
drag which the propulsor must overcome. By using consistent arrangement assump- 
tions and parametrics, the hull envelope sizes are reasonably correct relative to each 
other. Consequently, the propulsor characteristic provide a solid pairwise comparison 


between a full stern designed submarine and a tapered stern designed submarine. 
7.2.1 Need for Detailed Arrangements Study 


The parametrics level of arrangements consider both the submarine weight and vol- 
ume. While historically submarines have been weight limited, this study considers 


that future submarines may be volume limited as electronics become smaller and 
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LBP (ft) 
Sub. Disp. (ltons) 
SHP (hp) 
D (ft) 
L/D 
Volume (ft?) 106,800 
Wet. Surf. (ft?) 15,140 17,800 
Vol/ Wet Surf 6.91 
ER length (ft) oy 
ER volume (ft?) 17,080 
Propulsive Coefficient Ore" 
Estimated Max speed (knots) 21.8 
Stern Shape Full Stern Tapered 
Features PWR PWR 
Elect Dr Elect Dr 
Ext Motor Ext Motor 
* Note: This PC is a best estimate value which allows the speed 
of a full stern submarine to be estimated using traditional naval 
architecture scaling. It does not equal the actual PC, defined as 
PC = £45, due to separation near the stern when the full stern 


SHP’ 
body is towed. 


Vortex Lattice / RANS | Full Stern 


Body Press Drag Coefficient 
Duct Press Drag Coefficient 
Total Press Drag Coefficient 
Body Viscous Drag Coefficient 
Duct Viscous Drag Coefficient 
Total Viscous Drag Coefficient 
Stator Drag Coefficient 

Total Model Drag Coefficient 
Full Scale Correction 

Full Scale Drag Coefficient 
Advance Coefficient (Js) 
Thrust Coefficient (Ay) 
Torque Coefficient (Kg) 
Power Coefficient (Cp) 

Max speed (knots) 











Table 7.1: Summary of Variant 3 and 4. 
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lighter. The parametric study takes the most limiting of the two conditions while 
recognizing that, in the end, the weights will equal the volume. 

None the less, the need for a detailed study on area, volume and weight must be 
done to fully quantify the arrangement issues of a full stern design. Reference [11] 
does this on a more detailed level, but still at the preliminary design level. Until a 
design goes into the detailed design phase the real answer to the arrangements will 


remain unclear. 


7.3 Hydrodynamic Optimization 


7.3.1 Methodology 


The full stern variant could certainly be made better through more optimization. 
However, the quantity of variables that exist 1s considerable. The conclusion must 
be that the final Variant 3 hydrodynamic performance is at a local optimum. Strong 
interaction of the stern, duct, stator and rotor make the global optimization extremely 
difficult. 

On the other hand, the Variant 4 has relatively weak interaction of the stern, 
duct, stator and rotor. The stern taper allows the propulsor to be designed ignoring 
interactions, then superimpose the propulsor characteristics on the stern and iterate 
to reach the global optimum. 

A comparison of the hydrodynamic performance of the two variants yield insight 
on the full stern performance. Figure 7-1 shows the integrated drag due to pressure 
and viscous influences. The full stern has a large pressure drag near the entrance to 
the propulsor. However, the pressure drag is recovered as the flow enters duct. As 
discussed in Chapter 5, the body profile and its derivatives play an important role 
in the recovery of the pressure drag as well as the interaction between the duct and 


body. 
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Figure 7-1: For both the full stern (y = tanh(xr)) and ducted-tapered stern bodies, the integrated 
drag due to pressure and viscous influences are shown. 


7.3.2 Need for Hydrodynamic Tool for Full Stern Parametric Studies 


In optimizing the hydrodynamics performance of the full stern, only a limited number 
of cases are possible. The strong interaction of the stern components forbids a designer 
from optimizing one component and stating that the whole propulsor is optimized. 
Therefore, a design tool needs to be invented that can optimize the entire propulsor 
at once. 

Existing blade design methods for marine propulsors. such as PLL developed 


at MIT [6], are very often based upon preliminary lifting line procedures. These 





represent a propeller by a set of straight radial lifting lines and the wake as purely 
helical, i.e. vortices trail aft on constant radius cylinders. The limitation of these 
assumptions are reached since the full stern vehicles have geometries with a rapid 
diameter decrease at the after end. The straight radial lifting lines and helical wakes 
are inappropriate. 

One alternate preliminary design technique which models contracted stream tubes 
geometries is a streamline curvature method [22]. A shortcoming of streamline cur- 
vature methods, however, is that they typically do not model the flow outside of the 
duct. The duct streamline is used as the outer flow boundary. Therefore, the effect 
of the duct is not completely included. 

Work towards an improved propulsor design method is underway [27], [19], [28]. 
The new technique, termed DPLL for Ducted Propeller Lifting Line, is under devel- 
opment at MIT. The benefits should be significant as a typical case can be analyzed in 
one tenth of the time it takes now. The computational advantages include the ability 
to completely understand and model a full stern submarine to achieve the maximum 
performance. 

Work should be continued on these tools and others so that the maximum perfor- 
mance of a full stern body can be quantified. And, if appropriate, a full stern body 


should be designed and built to correlate theoretical predictions with actual results. 
7.3.3 Future Studies 


In reviewing the methods and the results in this thesis, several ideas present them- 
selves which may be of use in future studies. As discussed previously, the vortex 
lattice / RANS technique requires so much computational time that analyzing all 
possible cases is not practicable. By studying a carefully selected set of bodies, the 
powering trends can be identified in order to provide the basic data for a future full 


stern submarine. The following suggestions are made: 
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. Reduce Computational Time. Use techniques that calculate less detail at 
the expense of accuracy. Gauge the overall trends by using more cases with less 


accurate methods. 


. Full Verses Tapered Body Comparison. When comparing a full stern de- 
sign to a tapered stern design, the bodies must represent reasonable submarines. 
Comparisons based on equal volume are inappropriate since it may not be pos- 
sible to construct a submarine within the volume constraint due to stack length 


issues. 


. Quantitative Comparisons. Performance comparisons should be based on 
the maximum speed for a given installed shaft horsepower. A non-dimensional 
quantity which represents this comparison is the power coefficient. Compar- 
isons between a full stern and tapered stern using the propulsive coefficient are 


extremely misleading. 


. Global Optimization. Simultaneous optimization of the entire submarine, in- 
cluding the propulsor, is desired. However, in present day, it cannot be done. To 
design a full stern submarine, the ideas used in this thesis should be extended as 
discussed in Chapter 5. Start with a parent hull form, then meticulously examine 


daughter hull forms so that each daughter varies one parameter of interest. 


. Parent Hull Form. Chose a parent hull form which has the characteristics in 
the middle of the expected evaluated range. If the parent is in the middle of the 
range to be evaluated, then the trends developed from the daughters will likely 
hold. 


. Optimization Techniques. The parent-daughter method is analogous to a hill 
climber optimization technique. Other methods will work, however, no computer 
tool currently exist to vary the all of the propulsor characteristics which are of 


interest to the designer. Manual intervention is necessary to ensure a complete 
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10. 


EE 


sampling of characteristics. 


. Body Parameters. From forward to aft, the full stern submarine can be de- 


scribed by: nose radius, length forward (forward perpendicular to the maximum 
body diameter), length of parallel midbody, length aft (length from maximum di- 
ameter to the aft perpendicular), radius at shoulder, inflection radius and radius 


at tailcone. 


. Duct Parameters. The duct characteristics which may affect powering in- 


clude: thickness/chord, thickness distribution, camber/chord, camber distribu- 
tion, chord length, angle of attack, axial location relative to body and radius of 


inner surface at a reference location. 


. Propulsor Parameters. The parameters associated with the propulsor which 


may affect powering include: number of pre-swirl blades, axial location and 
characteristics, number of rotor blades, axial location and characteristics, and 


the tip clearance between the rotor and duct inner surface. 


Towed Verses Propelled. In a tapered afterbody, the body can be optimized 
for minimal drag. The propulsor is then added and the results are close to a 
global optimum due to the relatively weak interaction of the propulsor com- 
ponents. In a full stern, which will separate without an operating propulsor, 
optimizing the body without an operating propulsor would yield little relevant 
data. Instead, the optimization analysis must include an integrated propulsor to 


minimize powering. 


Body Continuity. A surprising result of Chapter 5 was the marked reduction 
in drag due to only slight afterbody changes. By designing the afterbody with 
smooth higher derivatives, the powering improved. To capitalize on this effect, 
the parent and all daughter bodies must have a smooth transition from the 


parallel mid-body to the propulsor. 
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13. 


Body Description. In order to maintain smooth bodies, it may be possible 
to represent the body profile as a power series or as a B-spline. The number of 


constraints placed on the body profile dictates the order of the curve. 


Computer Model To Physical Model. Lastly, testing of a physical model 
should be done to provide accurate and quantifiable results which can be used 
to validate computer models. Once, a full stern computer model is completely 
validated, the full stern series can be done to optimize a full stern body with 


great confidence. 
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